Anisotropic Dependence of Superconductivity on Uniaxial Pressure in Celrlns 
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We measure the effect of uniaxial pressure on the superconducting transition temperature Tc in Celrlns. We 
find a hnear change in Tc with both a-axis and c-axis pressure, with slopes of 56 mK/kbar and -66 mK/kbar, 
respectively. By comparing results from doping studies and different types of pressure measurements, we sepa- 
rate the influences of hybridization and dimensionality on Tc. We find the true geometric influence, for constant 
hybridization, is dTc/d(c/a) = 44 K. 



The low-temperature phases of heavy fermion materials 
have attracted much attention in recent years. They ex- 
hibit a range of correlated phases, including several types 
of magnetism. Superconducting regimes emerge near zero- 
temperature magnetic phase transitions. Non-Fermi liquid be- 
havior also appears near these quantum critical points, and can 
persist to significant temperatures. Tuning through alloying, 
pressure, or applied field allows exploration of the exact bal- 
ance among the phases. 

One of the most-studied heavy fermions recently has been 
CeMIn5 (M = Ir, Rh, Co). The proximity of these Ce- 
based 115 compounds to an antiferromagnetic quantum criti- 
cal point leads to a rich phase diagram CeRhlns at ambi- 
ent pressure has a superconducting transition near 0. 1 K, deep 
within an antiferromagnetic phase f3l. Hydrostatic pressure 
destroys the magnetism and raises Tc to a maximum of 2. 1 K 
at 16 kbar ^ 0Q. On the other hand, Celrlns and CeCoIns 
superconduct at ambient pressure iH [^, with several simi- 
larities to the high-temperature cuprate superconductors. The 
115 materials have a tetragonal, HoCoGas crystal structure 
which can be viewed as alternating layers of Celns and MIn2 
stacked along the (0 1) direction ol HI |^, reminiscent of 
the copper-oxygen planes in the cuprates. Power-law temper- 
ature dependences in heat capacity, thermal conductivity and 
spin-lattice relaxation rate in the superconducting state sug- 
gest c?-wave symmetry of the superconducting order param- 
eter 101. Furthermore, emergence of superconductivity near 
T/v and coexistence of homogeneous antiferromagnetism and 
superconductivity provide evidence that magnetic fluctuations 
mediate Cooper-pairing in these systems ISl llOlfTlll . 

One major influence on the superconducting transition tem- 
perature is the hybridization between the Ce /-electrons and 
In p-electrons lfl2i[l3i[l4ll . which controls the spin-fluctuation 
temperature, Tgf. Tsf is proportional to Tc and roughly in- 
versely proportional to 7, the Sommerfeld coefficient of the 
normal-phase specific heat. Smaller 7's in the isostructural 
PuMGas coumpounds (M = Rh, Co) partially account for their 
much higher superconducting transition temperatures. How- 
ever, even within the CeMIns series, 7 alone cannot describe 
all the variation in T^, for example, while 7 changes by just 
over a factor of 2 between M = Ir and M = Co, Tc changes 



by nearly a factor of 6. Thus hybridization cannot be the sole 
influence on Tc- 

Mean-field theoretical models of magnetically mediated su- 
perconductivity indicate a strong dependence of Tc on dimen- 
sionality ||ll,[lB[l3]. Measurements on CeMi^^jM^Ins show 
a linear relationship between Tc and the ratio of the tetrago- 
nal lattice constants c/a Ht]. Interestingly, for the Pu-based 
115 materials Tc is also linear in c/a, with the same relative 
change in Tc with dimensionality, ^ J^^[fa) ^^^^ ■ agree- 
ment in slopes makes sense if the difference in hybridization 
sets the overall temperature scale for each family but dimen- 
sionality governs the behavior within each family. Another 
way to control dimensionality and hybridization is by apply- 
ing pressure. Under hydrostatic pressure, c/a is not even 
monatonic for CeRhlns and CeCoIns, and Tc is not linear 
in c/a. However, all the CeMIns compounds have similar 
hybridization at the pressure Pmax which maximizes the Tc- 
Considering Tc and c/a at Pmax, where the effects of hy- 
bridization differences are reduced, does give a linear rela- 
tionship lfl3[ ITsll . Kumar et al. argue that hydrostatic pres- 
sure mainly alters the hybridization. The similar hybridiza- 
tions at Pmax suggest that hybridization determines the pres- 
sure Pmax while dimensionality governs the value of Tc at 

Pmax Cf. 

Uniaxial pressure is a natural technique for further explor- 
ing the effects of dimensionality, since the c/ a ratio of an indi- 
vidual sample can be increased or decreased depending on the 
pressure axis. Uniaxial pressure leads to fairly small changes 
in hybridization, since lattice constants decrease along the di- 
rection of applied force but increase in the perpendicular di- 
rections. On the other hand, the effects on c/a are much larger 
than for a similar hydrostatic pressure. To date, the effects of 
uniaxial pressure on the 115 materials have not been measured 
directly, although thermal expansion measurements combined 
with Ehrenfest relations predict the change in Tc with uniaxial 
pressure in the zero pressure limit lfl9ll . Here we explore the 
dependence of Tc on pressure and hence on c/a, and combine 
our results with those on hydrostatic pressure and alloying to 
extract the dependence of Tc on dimensionality. 

Samples were grown in alumina crucibles from a molten 
metal flux containing stoichiometric amounts of Ce and Ir, 



2 



and excess indium. We confirmed through transmission Laue 
X-ray diffraction that the lattice constants of our samples 
matched those reported elsewhere ll20ll . 

Our uniaxial pressure apparatus, shown in Figure [T] is a 
helium-activated bellows mounted on a dilution refrigerator 
ll2lll . and permits changes in pressure without thermally cy- 
cling the sample. We monitor the pressure through a piezo- 
electric crystal in the pressure column. The maximum achiev- 
able pressure depends on the size of the sample, but is typi- 
cally about 10 kbar. We measure the superconducting transi- 
tion with adiabatic heat capacity. Superconducting NbTi spac- 
ers between the sample and the pressure cell serve as the ther- 
mal link. 
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FIG. 1: Helium bellows setup for measuring heat capacity under uni- 
axial pressure. 



As can be seen from Figure [U the sample is unconstrained 
in directions perpendicular to applied pressure. To apply pres- 
sure along a specific axis, we first orient our sample using 
Laue X-ray diffraction. We then polish the sample according 
to the desired orientation. We preserve as much bulk as pos- 
sible during the polishing, since the observed time constant 
for thermal decay is proportional to the sample's mass. Our 
largest samples for pressure along the c and a axes were ap- 
proximately 30 mg and 80 mg, respectively. 

Figure |2] shows representative heat capacity data. Pressure 
applied along the c axis shifts the transition to lower temper- 
ature, while a-axis pressure has the opposite effect. The data 
are scaled, with a single scaling constant used for each sample, 
regardless of pressure. The average C/T over all pressures at 
800 mK is set to 700 mJ/mol K^. 

For a-axis pressure, we fit a function 7^ + aiT" + a2T~'^ to 
the heat capacity in the superconducting phase jst). The form 
fits well at all pressures, with the exponent 7i always close to 
1, as expected for line nodes in the energy gap. For c-axis 
pressure, the reduced Tc leaves too small a temperature range 
for reliable fits. If we fix n = 1 and refit the data, we find a 
slight decrease in both 7s and ai with pressure, as opposed to 
the sharp change in 7^ observed with uniaxial ll2lll and hydro- 
static ll22ll23ll in other heavy fermion superconductors. 

We define Tc as the temperature where C /T equals the av- 
erage of its maximum value in the superconducting region 



and its value in the normal phase at the onset of the transi- 
tion. We also model the transition as an abrupt discontinuity 
with the restriction that the normal-phase entropy is the same 
as for the actual data. The two methods give consistent val- 
ues of Tc- However, under c-axis pressure Tc decreases un- 
til eventually there is not enough data in the superconducting 
region to extrapolate the value of C/T in an equal-entropy 
calculation. For consistency between c and a-axis pressure, 
all values we report here use the average-C/T method for Tc 
in both pressure directions. The transition temperatures for 
different pressures, obtained using this average C/T method, 
are shown in Figure [3] In our pressure range, the change in 
Tc is linear It equals 56 mK/kbar for a-axis pressure, -66 
mK/kbar for c-axis pressure. These values agree fairly well 
with those derived for the zero-pressure limit from thermal 
expansion data: 54 mK/kbar and -89 mK/kbar, respectively 
III9I1 . For the tetragonal crystal structure, we can also com- 
pare our results to the effect of hydrostatic pressure through 
dTc/dpv = 2dTc/dpa + dTc/dpc- Our data yield a value 
of dTc/dpv = 46 mK/kbar, somewhat larger than the 25 
mK/kbar obtained through direct hydrostatic pressure mea- 
surements 1I24I1. 



We also comment on the normal state heat capacity. We 
find no significant change with a-axis pressure but an increase 
0.03 ± 0.01 J/kbar mol K? with c-axis pressure. Previous hy- 
drostatic pressure measurements ll24ll found a decrease of - 
0.02 J/kbar mol K^. Although the sign differences here and 
in STc/SP may be related, the magnitudes of the heat ca- 
pacity changes do not correspond to the overall Ce-In hy- 
bridization for the different types of pressure. Instead they 
suggest that hybridization in certain directions affects the heat 
capacity more than in others. The superconducting jump av- 
erages AC/'jTc = 0.73 ± 0.05 J/mol K^ and has no appar- 
ent trend with pressure, consistent with previous experiments 

We take as the width of the transition the span between 









p 


=0.62 kbar 






o-e> 


p 


=0.96 kbar 


1.2 






p 


=2.21 kbar 



o 1 
S 



0.8 



0.6 - 



•-•P 


=0.10 kbar 


□-qP 


= 1.69 kbar 




=3.71 kbar 




400 600 
Temperature (mK) 

FIG. 2: Heat capacity at various pressures. The data have been scaled 
by the normal phase heat capacity; see text. 
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the temperatures corresponding to 20% and 80% of the C/T 
range during the transition, as shown in the inset of Figure 
[3] Under pressure the transition becomes broader and more 
rounded. Ahhough an inhomogenous pressure distribution 
across the face of the sample or a variation in the cross- 
sectional area over the height of the sample would broaden 
the transition, reasonable values for these effects would pro- 
duce less than 15% of the observed broadening. 
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FIG. 4: Tc vs. c/a. For a-axis pressure, we calculate c/a using the 
minimum, geometric mean, and maximum a values; see text for de- 
tails. The solid line is a least-squares fit to the c-axis data, with slope 
35 K. The dashed line is a least-squares fit to the a-axis data using the 
mean values for a. Inset: putative Tc vs. c/a from substitution and 
pressure data, with hybridization potential Vp / equal for all points. 



FIG. 3: Transition temperature as a function of uniaxial pressure. Tc 
was determined using the temperature corresponding to the average 
of the C/T values at the peak and onset of the transition. Lines are 
linear fits, with slope 56 mK/kbar and -66 mK/kbar. Inset: 20%-80% 
transition width vs. pressure. The dotted line is a guide to the eye. 

We apply a small pressure during cooldown to keep the 
sample in place. This initial pressure differs for each 
cooldown but is typically on the order of 0.3 kbar. To deter- 
mine the initial pressure, we extrapolate Tc versus p to where 
the transition temperature reaches that measured on samples 
outside of the pressure cell. The pressure labels for the curves 
in Figure |2] are adjusted for the initial pressure, so that the 
values listed are the correct pressures. 

Figure |4]redisplays the data of Figure [3]to Tc as a function 
of c/a. We use room temperature lattice constants fl^ and the 
elastic constants of CeRhlns [18] to compute c/a at different 
pressures luOll. We assume that the Celrln.g elastic constants 
are similar ||19ll . Also, at our low to moderate pressures, the 
stress/strain relations are still linear. 

Just as for measurements across the CeMIn.g family, we find 
that Tc increases with increasing c/a. The dotted line of Fig- 
ure|4]is at the ambient-pressure value for c/a. Applying c-axis 
pressure decreases c while increasing a through the Poisson 
ratio, so the data from c-axis pressures appear to the left of the 
dotted line. On the other hand, a-axis pressure decreases the 
a lattice constant along the pressure direction and increases 
c, but also increases the a-axis lattice constant perpendicular 
to the pressure direction. Since the largest change in lattice 
constant is along the pressure direction, the overall effect is 
an increase in c/a, and the a-axis data appear to the right of 
the dotted line. However, for a-axis pressure the ratio c/a is 
not a single well-defined quantity. The value of a is minimum 



in the direction parallel to the applied pressure and maximum 
perpendicular to the pressure. Both of these extreme values 
are used to calculate c/a, with the results plotted in Figure]?] 
The remaining data set in Figure]?] uses the geometric mean 
of the maximum and minimum a values, which is the natural 
way to compare in-plane areas with the perpendicular lattice 
constant. As seen in Figure ]4j using the mean value gives 
a kink in dTc/d{c/a) at the crossover from c-axis to a-axis 
pressure. The larger slope of the a-axis data may indicate the 
influence of hybridization. Since the hybridization depends 
on the spacing between atoms, it increases for both directions 
of uniaxial pressure. However, the increase is small because 
the decrease in atomic spacing along the direction of the ap- 
plied pressure is partly compensated by increased spacings in 
the perpendicular directions. 

By considering several different methods of shifting Tc, 
we can separate the true influence of geometry on the tran- 
sition temperature from the effect of hybridization. To do 
this, in Table ]l] we consider the effects of hydrostatic pres- 
sure, chemical substitution, and our a-axis and c-axis pres- 
sure measurements. We use Harrison's calculation fZSl of 

Vpf in a tight binding approximation, ijpf ^ -^^-g?-^ ■ Here 
ijpf ~ 10-\/2T/ TT is a dimensionless constant, m is the electron 
mass, and rf = 0.445A and rp = 19.1 A are wave function 
radii llTsi 12611 . For rp we extrapolate values from ll26ll . p. 644. 
For hydrostatic pressure data, we take the percentage changes 
in Vpf and c/a from Table II and Figure 2 of lfl3ll . 

Each experiment in Table]!] — substitution, applied pressure, 
or a combination of the two — increases Vp / over its value in 
Celrln5 at ambient pressure, with a percent change given by 
AVpf. For each type of measurement we simultaneously scale 
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TABLE I: Values used in adjusting for effect of hybridization changes on Tc. Zero-pressure structural parameters from 12(1 12711 : hydrostatic 
pressure data from 1 13]. See discussion in text. 
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the changes in Tc, c/a, and Vpf hy a single factor to reach 
AVpf = 0.2%. Our scaled values show how Tc would vary as 
a function of c/a with Vp / held constant. The results from the 
final two columns of Table |T] also appear in the inset of Figure 
|4] The solid line is a best fit to all six points, with slope 44 
K. Significantly, the slope remains 44 K if the fit omits the c- 
axis pressure data, which has particularly small hybridization 
change and lies far from the other points. 

The broken symmetry perpendicular to the c-axis when 
pressure is applied along the a axis has little effect on the su- 
perconductivity, judging from the linearity of the fit in the Fig- 
ure |4]inset. We note that using the minimum value of the lat- 
tice constant a preserves the linearity of Tc vs. c/a. Although 
probably coincidental, this may indicate how the destruction 
of tetragonal symmetry affects the superconductivity. 



Our measurements show that up to a few kbar Tc changes 
linearly with uniaxial pressure, at 56 mK/kbar for a-axis pres- 
sure and -66 mK/kbar for c-axis pressure. This confirms the 
important role of c/a in controlling the onset of superconduc- 
tivity. Changes in fp hybridization between the Ce and neigh- 
boring In atoms also have a large effect. By comparing several 
techniques that alter the dimensionality and the hybridization 
to different degrees, we control for hybridization changes and 
find a pure geometric influence of dTc/d{c/a) ~ 44 K. 
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